Fermilab-Pub-09-662-T 



Spin Correlation Effects in 
Top Quark Pair Production at the LHC 

Gregory MahloEB 

Penn State Mont Alto 
1 Campus Drive, 
Mont Alto, PA 11231, USA 

Q ■ Stephen J. Parkfl 

(N 

I Theoretical Physics Department 

Fermi National Accelerator Laboratory 
P.O. Box 500, Batavia, IL 60510, USA 



(N 



in 

(N 

, (Dated: January 17, 2010) 

'p^' Abstract 

1^ , At a 14 TeV proton-proton collider, the Large Hadron Collider (LHC), we show that top quark 

pair production is dominated at low invariant mass by the fusion of two like-helicity gluons, pro- 
ducing top quark pairs in the left-left or right-right helicity configurations. Whereas, at higher 
^ ' invariant mass the production is dominated by the fusion of unlike-helicity gluons, producing top 

■ quark pairs in the up-down or down- up off-diagonal configurations, identical to top quark pair pro- 

T^lj- . duction via quark-antiquark annihilation. We study in detail the low invariant mass region, and 

\ show that the spin correlations can be easily observed in this region by looking at the distribution 

of the difference in the azimuthal angles, A(/), of the dileptons decay products of the top quarks in 
f-*) ' the laboratory frame. Due to the large cross section for top pair production at the LHC, even with 

a cut requiring that the invariant mass of the top quark pair be less than 400 GeV, the approx- 
imate yield would be 10^ di-lepton (e, //) events per fb~^ before detector efficiencies are applied. 
^ , Therefore, there is ample statistics to form the distribution of the dilepton events, even with 

^ \ the invariant mass restriction. We also discuss possibilities for observing these spin correlations in 

the lepton plus jets channel. 

PACS numbers: 14.65.Ha 



* Electronic address: 'gdmlO@psu.edu 
^Electronic address: parke@fnal.gov 



1 



I. INTRODUCTION 



Since the discovery of the top quark with a mass between 150 and 200 GeV^ by the CDF 



and DO experiments at Fermilab in 1995, numerous authors, see [2[-[10|, have asked the 
question "Can the spin correlations in top quark pair production be observed?" This is a 
vahd question since the top quark hfetime in the Standard Model is very short compared to 
the spin de-correlation time for such a heavy quark. In particular, Fj- ~ Gpfn^ » AQQ^/rrit 
where Ft is the total width of the top quark, rrit is mass of top quark, Gp is the Fermi 
constant and Aqcd is the QCD scale, thus the top quark decays before QCD interactions 
have the opportunity to appreciably affect its spin. The angular distribution of the top 
quark decay products in t — )■ + b followed by + u or d + u are correlated with 

the top spin axis as follows: 

I ij^Y f +1.0 or d-quark 

1^1 + ctj cos Xi)/2 Cii = \ —0.31 z/ or u-quark (1) 



^Td cos Xi 1 _o.4i 5.quark 



where Xi is the angle between the i-th decay product and the top quark spin axis in the 
top quark rest frame. Clearly, the charged lepton or (i-quark coming from the decay of the 
VT-boson are the most correlated with the top quark spin axis. For the anti-top, the signs 
of the ttj coefficients are flipped. Thus, if the spins of the top are correlated in top quark 
pair production and since the decay products of the tops are correlated with the spins then 
the decay products of the two top quarks are correlated. Since there is no net polarization 
of the top quarks, at least to leading order, the correlation between the i-th decay product 
of the top and z-th decay product of the anti-top can be expressed by 



1 rfV _ 1^ 

ax d cos Xid cos Xi 4 

with 



1 + Ctt ai aj cos Xi cos Xi) ■ (2) 



cx^/l tA i^ production cross section for top quark pairs where the top quark has spin up 
or down with respect to the top spin axis and the anti-top has spin up or down with respect 
to the antitop spin axis. Clearly, the right choice for the spin axes of the top quark pair 
is important since a poor choice of spin axes can lead to a small value for the correlation 
parameter, Cu and hence to small correlations between the decay products of the top and 
the antitop. For some processes, e.g. qq — )■ ti, there exist spin axis choices such that the 
correlation parameter, Cu, is maximal. However, even with these optimal choices, to observe 
the correlations one has to measure the angles Xi and Xi between the decay products and 
the spin axes in the rest frame of the top and antitop quarks respectively. To do this one 
has to reconstruct the top and the antitop rest frames; this is very challenging at a hadron 
collider. An obvious question is "Are there variables that carry the signature of the spin 
correlations which can be measured in the laboratory frame?" For qq — )■ tt no such variables 



^ The current Tevatron-averaged, best-fit value for the top quark mass is 173.1 ± 1.3 GeV, see Ref. Jj. 
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have been found that show a significant difference between full spin correlations and no spin 
correlations. 

It has been known for some time that the spin correlations at the LHC are described 
very well by the helicity basis at low top-antitop invariant mass, whereas at higher invariant 
mass, counter to one's naive expectation, the spin correlations are degraded in the helicity 
basis. In this paper we explain this phenomenon by showing that at low invariant mass top 
quark pair production is dominated by like-helicity gluons, producing top quark pairs in the 
left-left or right-right helicity configuration, independent of the invariant mass. Whereas, 
at higher invariant mass top quark pair production is dominated by unlike-helicity gluons, 
producing top quark pairs in the up-down or down-up off-diagonal configuration, identical to 
that of top quark pair production via quark-antiquark annihilation. At ultra-high invariant 
masses, the up-down and down-up off-diagonal configurations become the familiar left-right 
or right-left helicity configurations; however, at the LHC only a small fraction of the total 
number of produced top pair events are in this ultra-high invariant mass region. The fact 
that the contributions from like and unlike helicity gluons impart different spin correlations 
to the top quark pairs makes gg — )■ tt, which dominates top quark pair production at the 
LHC, a much richer process for analysis than the dominant Tevatron mechanism qq — )■ ti. 

We study the low invariant mass region in detail. In this region, the like-helicity gluons 
dominate the production. We show that the spin correlations can be easily observed in 
this region by looking at the distribution of the difference in the azimuthal angles, A0, of 
the dilepton decay products of the top quarks. For top quark pairs with an invariant mass 
< 400 GeV, the spin correlations give a 40% enhancement of this distribution at small angles 
(A0 !=s 0) and a 40% suppression of this distribution at large angles (A0 ^ tt) compared to 
if there were no correlations between the production and decay of the top quarks. About 
20% (~ 200 pb) of the total next-to-leading order top-antitop quark production cross section 
passes this invariant mass cut so even in the di-lepton channel there are large numbers of 
events available to measure this distribution: about 10^ events per fb~^ before detector 
efficiencies are applied. 

The outline of this paper is as follows: in Sect. [Ill we review what is known about the spin 
correlations for qq — )■ ti before taking a closer look at the spin correlations for gg — >■ tt in 
Section IIIII Section [IV] addresses the question of how the top quark pair events at the LHC 
populate the scattering angle versus invariant mass plane. In Section |V] we add the decays of 
the top quarks. In Sect. I VII we address the issue of which angular distributions are sensitive 
to the presence or absence of angular correlations among the tt decay products. We show that 
at low invariant mass the difference in the azimuthal angle of the charged leptons carries 
the signature of spin correlations in the laboratory frame. We also discuss possibilities 
for observing these spin correlations in the lepton plus jets channel. We summarize our 
conclusions in Sec. IVIII Lastly, we include an Appendix which outlines the highly efficient 
method used in this paper for calculating the spin amphtudes. It can be applied to any 
2 — 2 process with arbitrary spins of the final state particles. 
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II. REVIEW OF qq tt 

A. Spin Amplitudes in an Arbitrary Basis 

For a massive particle with momentum, t, and spin vector, Sj, the following relationships 
are satisfied 

= m^, St = -1, and t-St = 0. (4) 

The spin vector st is most conveniently defined in the rest frame of the massive particle; in 
this frame it only can have spatial components. Thus, for top-antitop quark pair production 
via quark- ant iquark annihilation or gluon-gluon fusion, we define the spin vector st in the 
rest frame of the top quark. Since CP is conserved at tree level for this process, we restrict 
the spin vector to be in the scattering plane. It is convenient to measure the direction of 
spin vector st with respect to the antitop quark or recoil direction. Thus we define the unit 
vector St such that its direction is at an angle ^, measured clockwise with respect to the 
antitop quark direction (see Fig. [T]). 

For the antitop quark we proceed in a similar fashion in defining St- Here, instead of 
using the same angle ^ to specify the direction of the spin vector with respect to the recoil 
direction, we use a different angle This allows for the independent manipulation of the 
antitop and top quark spins at intermediate steps. However, in the end, we will set ^' = ^: 
then, in the zero momentum frame (ZMF), the spin vectors of the top and antitop quark 
are back-to-back as expected from the symmetry arguments. 

It is particularly convenient to use the following combinations of t and St as follows: 

ti = {t + mtSt)/2, and t2 = {t - rritSt) /2. (5) 

These satisfy 

tl = tl = {], t = ti+ 12, and 2ti ■ t2 = m]. (6) 

{ti and ^2 are defined similarly for the antitop quark.) Since ti and t2 are light-like vectors, 
the full power of the spinor helicity method can be used in evaluating the amplitudes as 
discussed in Ref. fnl, resulting in many simplifications. For example, the Dirac spinor for 



a massive fermion with spin up can be written 



U,{t) = ^f/(ti) + e^^l^Uit2). (7) 

This factorization into chiral components, with one depending only on ti and other only 
t2 (apart from the phase factor e** = U {t2)^^^U {ti) / rrit) , is particularly useful. For our 
purposes in this paper we need only two such spinors: a spinor for the top quark, and a 
spinor for the antitop quark. These are given in spinor helicity notation by 

Ht) = (ti + I + ^ - I = — (ti + |(t + m,) 

nit 

Vtii) = - \i2+)^-^^^^ = {mt - — (8) 

nit T^t 



The full set of spinors for massive fermions may be found in Ref. [2|. In this paper we will 
show that once you have the amplitudes for gg — )■ t^i^ in terms of the spin angles ^ and 
you can obtain all of the other spin combinations for the top and antitop by simple algebraic 
manipulations. 
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Top quark rest frame 




(recoil direction) I 

FIG. 1: Spin vector for top quark pair production. The direction of the top quark spin vector sj 
is given by the angle ^ in this frame, measured in the clockwise direction from the recoil direction. 
The inset illustrates the situation in the ZMF, where the top and antitop are produced back-to-back 
and the scattering angle 9 from the incoming beam direction. 

B. Review of the Spin Structure of qq tt 

The spinor structure of the matrix element for gnqL — > tt can easily be shown to be of 
the form 

Uit){\q+){q+\ + \q-){q-\}Vit) (9) 

by using the Fierz identities on the current- current structure of the matrix element given by 
the standard Feynman rules. This can be used to recover the well-known tree level matrix 
element squared for gnqL — ?• ti (see Ref. [1]): 

l-AiqnqL ^ t^i^ and = \A{qLqR t^i^ and t|t^)|^ 

~ (7^^ sin cos^ — cos ^ sin^)^ 

l^iquqi^t^ii or tii^)]"^ = \A{qLqR-^tit^ or t^t;)^ 

~ (7"-^ sin6'sin^ + cos6'cos^ =F 1)^. (10) 

Here the same spin angle, ^, has been used for both the top and antitop quarks. From this 
general result it is clear that there is a basis, defined by 

tan^ = 7"4an6', (11) 

which sets the + ii component to identically zero for all /3, leaving only the tt + tt 
component. This basis was first identified by Parke and Shadmi in Ref. [sj and has been 
called the off-diagonal basis. It interpolates between the beamline basis (cos^ = cos^) at 
threshold, and the helicity basis (cos^ = ±1) in the ultra-relativistic limit. In the off- 
diagonal basis, Eq. ( ITU]) becomes 

\-A{qRqL t^i^ and t^t;)!^ = \A{qLqR t^t^ and t^t;)]^ = 

lAiqnqL ^ t^i^ or t^t^)\^ = \A{qLqR ^ t^i^ or t^i^)\^ ~ (^1 T ^1 - /^^ sin^ 6^ , 

(12) 
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whereas the hehcity basis is obtained by setting cos^ = —1; in this basis 

\A{qRqL tntR and tLh^'^ = \A{qLqR tntR and tLti)]^ ~ 7"^sin^6' 

\AiqRqL tRh or titR)]^ = \A{qLqR ^ titR or tRh)]^ ~ (1 ± cos 0)1 (13) 

Clearly, for 7 >> 1, the helicity basis and the Off-Diagonal basis become identical. As we 
will see in the next section, the spin correlations for unlike-helicity gluons producing top 
quark pairs are identical to those in quark-antiquark annihilation. 




FIG. 2: The spin configurations for the process qRqi tt are best described by the off-diagonal 
basis which interpolates between the beamline basis at low /3 to helicity at very high /3 as given 
by Eq. [T5J (a) is the limit /3 — t- where the top quark spins are aligned in the same direction as 
the incoming quark spins whereas (b) is the limit /3 — )• 1 where the helicity state tRtL dominates 
for scattering angles less than 90 degrees. The relative probability of tRti to titR is given by 
(1 + cos 0)2 : (1 - cos 0)2. 



III. A CLOSER LOOK AT gg 



tt 



The tree-level matrix element for gg — )■ tt can be factorized into two terms: one depending 
on the color factors and t and w-channel propagators and the other depending on the the 
spin of the gluons and top quarks, as follows 



Aigig2 ti) 



(2t-pi) 



+ 



'111 



(2t-P2) 



M{gig2 ^ tt). 



(14) 



The reduced matrix element M{gig2 — ?■ tt) is symmetric under the interchange of the two 
gluon momenta but depends on the the helicity of the gluons and the spin of the top and 
antitop quarks. 

The square of the color-propogator factor, summed over the gluon and top quark colors, 
is given by 



E 

color 



{2t ■ pi) {2t ■ P2) 



4 4{t ■ p,)^ + A{t ■ p2)^ - {t ■ pi){t ■ P2) 
3 (t-pi)2(t-p2)2 



(15) 
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When evaluated in the ZMF, this sum reduces to the form 



color 

with 



{2t ■ pi) {2t ■ P2) 



^ ^ (16) 



W,c.) = j7T^^- (17) 



3 (1 - P'cj 

In these expressions (3 is the ZMF speed of the top quarks and cg is the cosine of the ZMF 
scattering angle 9. 

The reduced matrix element for on-mass-shell top quarks, M{gig2 tt), is simply given 

by 

MigngL ^ ti) = lSPl±}lh±lu{t){\p,+){p2 + | + \p2-){pi - \}V{i) (18) 

2pi ■ P2 

for unlike helicity gluons and by 

M{gngR^ti) = 2mt^^^^p^U{t)jLV{t) where 7^,^ = 1(1^75) (19) 

{Pi + \P2-) 2 

for like-helicity gluons. Note the similarity in the spinor structure for gngL — ?• tt and 
QrQl tt. Also, the spinor structure for gngn — )■ ti is particularly simple, U{t)'yLV{t); it 
contains no s-channel pole. In a later section of this paper, we use these two expressions 
to give a simple analytic expression for gg — )• ti including the decay of the two top quarks. 
However, in the next section we will evaluate these expressions using the spinors for polarized 
top quarks given in Eq. ([8]). 



A. Unlike-Helicity Gluons 

For unlike-helicity gluons the reduced matrix element M{gB,gL ^t^t) given by 

M{gRgL ^ tt^t) = + l^2bi+)(P2 + \ii-)-{ti + \p2-){pi - Mh-)], 

(20) 

which, when evaluated in the ZMF using the spin vectors described in the previous section, 
becomes 

M{gRgL t^i^) ~ /3sin^{(l - cos^) sin(^/2) cos(^V2) - (1 + cos^) cos(^/2) sin(^'/2) 

+7-' sin^[cos(e/2) cos(e72) - sin(e/2) sin(e72)]}. (21) 

Here the coefficients in front of the products of the ^-dependent trigonometric functions 
are the appropriate helicity amplitudes whereas the products of the ^-dependent trigonomet- 
ric functions themselves are products of Wigner d-functions (see Appendix A). The relative 
signs between the various components of these expressions are important and care must be 
taken to make sure they are correct. 
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Using a different spin angle for the t and t allows for manipulation of the spin of the top 
independent of the antitop and vice versa. Thus, all of the spin amplitudes for gngi — ^ ti 
can be simply obtained from QrQl — ^ ^t^t follows: 



|M^t(^,OI = |Mtt(e±vr,ni = 
|Mt;(e,OI = |Mn(e,e'±vr)| = 
|M^(e,Ol = |i^n(e±vr,r±vr)| 



d 



d 

dm 



d 



d 



(22) 
(23) 

(24) 



dil2) \di'/2^ 

Flipping the spin of a particle is accomplished by one of two equivalent methods: 

• Addition or subtraction of tt from the spin angle ^. 

• Differentiation of the amplitude with respect to ^/2. 

A detailed discussion with examples of how to use these techniques for arbitrary spins is 
given in Appendix A. 

At this stage we can make the spin axes of the top quark pair back-to-back in the ZMF 
by setting ^ = ^. Thus, for unlike-helicity gluons we obtain 



\M9r9l MM and 
\~A{gRgL -> t^t^ or t^i^ 



3^(/3, 6) 13' sin' ^(7^^ sin ^ cos ^ - cos 6 sin 



(25) 



y{/3, 6) j3 sin ^(7 sin ^sin^ + cos^cos^T 1)^ (26) 



and 



with 



\-A{9LgR^ and t^t^)^ = \A{gRgL ^ and tit^)f 
\A{gLgR^ht^ or t|f^)p = \A{gRgL^t^k hh)\'^^ 



\A9R9L ^ tt)\^ = \A9LgR ^ ti)\^ = 2 y{f3, 6) f3' sin^ 9{2 - f3^ sin' 



(27) 
(28) 



(29) 



all 



all 



As in qq — > tt, a great simplification occurs for the off-diagonal basis tan^ = 7 ^ tan^, 
where 



\AigLgR MM and t^t^ 



\A{gRgL t^t^ and t^t^)\' = 0. 



(30) 



The off-diagonal basis is the basis that interpolates from the beamline basis at threshold 
to the helicity bases at ultra-relativistic energies for the qq — t- tt process. Thus, for unlike- 
helicity gluons we have a very similar situation to that of qq — )■ tt: the only non-zero 
amplitudes are given by 



\A{gRgL Mh t^ttT = l^iOLgR ht^ or t^t^l 



y{/3,e) rsin'0 1^ -y/ 1-/32 sin' ^ 



(31) 



as illustrated in Fig. |3l 
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FIG. 3: The spin configurations for the process gngi — ^ tt are best described by the ofF-diagonal 
basis, which interpolates between the beamhne basis at low /3 and the helicity basis at very high ^ 
(see Eq. [3T]) . As far as the spins of the top quarks are concerned, this process, gnQL — ^ tt, is identical 
to top quark production via quark- antiquark collisions, gnQL — ^ tt, see Fig. [2j (a) illustrates the 
limit /3 — )• where the top quark spins are aligned in the same direction as the incoming gluon spins 
whereas (b) illustrates the limit /3 — >• 1 where the helicity state tRti dominates for scattering angles 
less than 90 degrees. The relative probability of tjitL to tLtR is given by (1 + cos 9)'^ : (1 — cos 9)'^. 

B. Like-Helicity Gluons 

For like-helicity gluons the reduced matrix element M{gRgii — )■ t-i-tf) is simply given by 
the following combination of spinor products 

Migngn -> t^t^) = 2m, ~ |^^+| (ti + (32) 

(Pl + \P2-) 

which when evaluated in the ZMF using the spin vectors described in the previous section 
is just 

M^gngn ^ t^t^) ~ 7^'{(1 - /?) cos(e/2) cos(e72) + (1 + /3) sin(e/2) sin(e72)}. (33) 

Treating these expressions in a manner similar to the unlike-helicity case discussed in the 
previous section we obtain 

\A{gRgR ^ t^t^ or t^tjp = y{l3,e) 7-^(1 ^ /3 cos 0', (34) 

\A{gRgR ^ t^i^ and t^^)\^ = y{P, 6) 7-^/32 gin^ ^. (35) 

Similarly, it is easy to show that for left-handed like-helicity gluons 

\A{gLgL -> hh or tt^t)r = \A9r9r ^t^t t^t^)\^ (36) 
\A{gLgL^hh and t^^)\^ = \A{gRgR^t^t^ and t^t^)\^ (37) 

Summing over all of the final spins gives 

l-^^angR -> ti)\' = J2 lAigigi ^ tt)\' = 2 y{f3,e) (1 - f3'), (38) 

all all 
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FIG. 4: The spin configurations for the process gngn — )• tt are best described by the hehcity basis 
for all /3 (see Eg. H0|) . (a) illustrates the limit /3 — )• where the relative probability of tRtu to t^tL 
is (1 + /3)^ : (1 — /3)^ whereas (b) illustrates the limit /? — )• 1 where tjitp; completely dominates. For 
9l9l ^ it, flip the spins on both the gluons and the top quarks. 



independent of the spin axis used for the top quarks. 

Clearly, a great simplification occurs for like-helicity gluons if one uses the helicity basis 
(^=0 or tt) for the top quarks. In the helicity basis 

\A{gLgL ^ tntL and t^tij)!^ = \A{gRgR tntL and titij)]^ = (39) 

for all values of f3. Conventional wisdom states that helicity provides a simple description 
for most processes only at ultra-relativistic energies. However, as illustrated in Fig. HJ ti 
production from like-helicity gluons is an exception to this expectation: in this case, the 
helicity basis provieds a simple description for all (3, with the only non-zero amplitudes 
given by 

\A{gRgR ^ tRtR or tLh)]^ = \A{gLgL -> or M^l' = 3^(/3, 0) 7"'(1 ± (40) 
Both the like and unlike gluon helicity amplitudes agree with those found in the appendix 



ofRef. 13 



C. Combining Like- and Unlike-Helicity Gluons 

At the LHC we must combine the like-helicity and unlike-helicity gluon cases since there 
is no way to polarize the incoming gluons. By looking at Eqns. ( 12 6 p and f l35p . it is clear 
that there is no basis which makes the top quark spins purely tt + ii OR ti + it at the 
LHC because the constant term appears in ft + H for like-helicity gluons and ti + it for 
unlike-helicity gluons. 

However, there are regions of the (cos 6*, /3) plane for which one of the like-helicity gluon 
amplitude or the unlike-helicity amplitude dominates. Along the curve given by 

sm6 = 1 or, equivalently, = 1/(2 - cos^ 6) (41) 

the like-helicity and the unlike-helicity contribute equally to top quark pair production. On 
this curve 

\AigRgR ^ tt)\^ = l-^idLgi ^ ti)\' = ^ \A{gRgL ^ tt)\^ = ^ lAigign -> «)!'• 

all all all all 
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In the region /37sin0 < 1 the hke-hehcity gluon amphtudes dominate the cross section, 
whereas in the region P'-fsinO > 1 the unhke-hehcity gluon amphtudes dominate the cross 
section. Thus, it is clear that one should use the helicity basis when P'fsinO << 1 and the 
off-diagonal basis when P'ysinO >> 1. In the next section we will optimize the basis choice 
to maximize the spin correlations in the intermediate region, /S'jsmO ~ 1. 



D. Optimizing the Choice of Spin Basis 

For unpolarized gluons, the fraction of top quark pair events at a given point in the 
(cos 9, /3) plane that have '|''|' or H spins is 



f{0,/3) ^ 



7-2(1 + /J^cos^^) + /32sin2 6'(7-^sin6'cos^ - cos6'sinO^ 
((1 - + /32 sin^ 0{2 - /32 sin^ 6)) 



It is a straightforward analytic exercise^ to find the extrema of this function with respect to 
the angle ^. The maxima, /same (^5 Z?), gives the maximum fraction of tt + ii whereas the 
minima, /oppo(6', /?), gives the minimum fraction of + li or, equivalently, the maximum 
fraction of + it- These fractions are given by 



/{same, oppo}(^)/3) — 

7-2 + insin^ 9 cos^ e + 7-2 sin^ 6 + T^) { 1 ± ^1 - } 

((1 - + /32 sin^ ^(2 - /32 sin^ ^)) 

Both extrema occur when ^ satisfies 



(42) 



27-1 sin^6' cos 6^ 

tan2(^|same, oppo} = . 2 n 9~2) _9 ■ 4 /i Zo' V^^) 

sm COS"' 6 — ^1 sm — 7 ^ 

they are related as follows: .^oppo = ■Csame + 7r/2. The contours of /same(6', /3) in the (cos 6, (3'^) 
plane are given by the solid lines in Fig. |5]^a) whereas for foppo{0,/3) see Fig. Mh)- 

At any given point in the (cos 6, P^) plane, the basis which exhibits the strongest correla- 
tions is the one whose spin fraction has the largest difference from i. If | /oppo (6', (3) — l/2\ is 
larger than |/same(6', /3) — 1/2| then one should use .^oppo; otherwise, .^same should be used. The 
condition that must be satisfied for both fsame{0,(3) and /oppo(6',/3) to have equal difference 
(but opposite sign) from 1/2 occurs when 

/same(^, P) + /oppo(^, P) = I, OT Pj siu 6 = 1. (44) 

Not surprisingly this is the same curve that also separates the dominance of the contribution 
of like-helicity from unlike-helicity gluons. Thus, when P'y sin 6 < 1 the like-helicity gluons 
dominate and ^same should be used to maximize the + ii fraction, whereas if P'j sin 6^ > 1 
the unlike-helicity gluons dominate and we should use ^oppo to maximize the + it fraction; 



^ The numerical solution was studied in Ref. 
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-1.0 -0.5 0.0 0.5 1.0 -1.0 -0.5 0.0 0.5 1.0 

COS 6 COS 6 



FIG. 5: Panels (a) and (b) show, in the (cos0,/3^) plane, the maximum fractions of ft + ii) 
/same(^, and ti + iti /oppo(^,/3), respectively (solid contours). In both panels the long-dashed 
line is /?7 sin 6 = 1 oi equivalently = 1/(2 — cos^ 6) below which the like-helicity gluons dominate 
and above the unlike-helicity gluons dominate. In (a) the short-dashed lines are the fractions of 
LL+RR top quark pairs in the helicity basis which below /3'y sin 6 = 1 are close to the maximum 
fractions for ft + /same(^, Whereas in (b) the short-dashed lines are the fractions of ti + it 
in the Off-Diagonal basis which above f3j sin 9 = 1 are close to the maximum fractions of tJ, + 

/oppo(^,/5). 

this is equivalent to minimizing the tt + ii fraction. The long-dashed line in both parts of 
Fig. |5]is the curve P'-fsinO = 1; this is the demarkation curve between maximizing tt + ii 
and maximizing fj, + if. 

It is worthwhile asking the following two questions: 

• In the region dominated by like-helicity gluons {P'jsinO < 1), how much does the 
maximum fraction of + il differ from what the helicity basis would give in the 
same region? 

• In the region dominated by unlike-helicity gluons [P'jsmO > 1), how much does the 
maximum fraction of + it differ from what the off-diagonal basis would give in the 
same region? 

These two questions are also addressed by Fig. [51 In (a) we have also plotted (short dashes) 
the fraction of top quark pairs which are LL+RR in the helicity basis and in (b) the fraction 
that are ti + it off-diagonal basis. Clearly, these figures indicate that the helicity 

basis does almost as well as the basis which maximizes tt + ii fo^' f^l sin ^ < 1 and that the 
off-diagonal basis does almost as well as the basis which maximizes ti + it fo^' P'ysinO > 1. 

Now that we understand the production dynamics for top quark pair production from 
gluon-gluon fusion we can turn to the question of what regions in the (cos 6',/?^) plane do 
the top quark pair events occur at the LHC. 
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TABLE I: Fraction of tt events catagorized by the helicities of the initial state gluons and location 
in the [cosO,/]'^) plane at the LHC with ^ = 14 TeV. 



P^S0 <1 /Sjse >1 all 



gg Like 55% 
gg Unhke 20% 
Total 75% 



10% 65% 
15% 35% 
25% 100% 



IV. PHENOMENOLGY OF TOP QUARK PAIR PRODUCTION AT THE LHC 

(l/g) d"a/(d^" dcg) 



0.8 



0.6 




0.4 



0.2 



0.0 



FIG. 6: The differential cross section (l/cJtot) d?cF/d{cos6) df3'^ for top quark pair production 
at the LHC assuming total beam energy of 14 TeV. The long-dashed line is /3'jsm9 = 1 (or 
= 1/(2 — cos^^)) is the demarkation line for the differential cross section to be dominated by 
like-helicity gluons (below) and unlike-helicity gluons (above). 



In Fig. |6]we plot the differential cross section 

1 rfV 



in the (cos 6',/?^) plane where Utot is the total cross section. This figure gives us the relative 
distribution of top quark pair events in this plane. A breakdown of the fraction of events 
for like- and unlike-helicity gluons broken into the appropriate regions in fi''ysQ are given in 
Table E 

Fig. [7] shows the differential cross section with respect to ICsamel, where 



C'same = 2/|-|-+4,^ — 1 and Coppo = 2/|;+4,| — 1 



(45) 
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2.0 



/' all 




same 



FIG. 7: The differential distributions of |Csame|i (1/o't) da / d\Cse,mc\- The dashed curves are for 
the optimal basis, Eqn. (j43p . whereas the solid curves employ the helicity basis when P^sq < 1 
and the off-diagonal basis when (3^sq > 1. The mean values of |Csame| for the different bases and 
regions are collected in Table HIl 

are the quantities which control the size of the correlations for any given spin basis (note 
that Coppo = — C'same)- In this figure we have split up the contributions into two pieces: one 
for /37sin0 < 1 and the other for /37sin^ > 1 . For /37sin0 < 1, we show the contribution 
in the basis which maximizes the ft + ]^\. component as well as the helicity basis; whereas 
for /37sin^ > 1, we show the contribution in the basis which maximizes t"i + V\ as well 
as the off-diagonal basis. This figure clearly shows that there are only small differences 
between using the best basis and the helicity basis for (3'-)sm6 < 1 and the best basis and 
the off-diagonal basis for P'-fsinO > 1. 

At the LHC, the total top quark pair production cross section is ~1 nb (at next-to-leading 
order) giving approximately 10^ ti per fb~^; therefore, significant cuts can be made on the 
data before the statistical uncertainties become comparable to the systematic uncertainties. 
Thus, we will concentrate on the low /3 region for the rest of this paper since in this region 
like-helicity gluons dominate the production cross section and the boost of the top quarks 
does not mask the spin correlations. 
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TABLE II: Mean values of |Csame| using either the optimal or appropriate choice of off-diagonal or 
helicity basis for events with different values of fi^sg at the LHC with ^/s = 14 TeV. 



Basis 






all 


Optimal basis 


0.55 


0.59 


0.43 


Helicity or off-diagonal 


0.53 


0.57 


0.39 



V. ADDING TOP QUARK DECAYS 



From Eqs. ( ITS]) and (fT^ . it is easy to add the decays of the on mass shell top quarks, 
t ^ h + e + V and i ^ b + fi + u, via the following replacements: 



9l 



V{i) {-i+mt)\fi+){iy+\b-) 



{b- \i^+){e+ \{t + mt) 



9l 



(46) 



The Fierz identity has been employed in the derivation of these replacements. Thus, the 
total matrix element squared for top quark production and decay via gluon fusion, summed 
over the colors of the incoming gluons and outgoing 6-quarks, is given by 



\A\l, + \A\ln = }C 



{2{pi ■ t){p2 ■ t) -mf{pi-p2) 



2(t-e)(t-/i) {pi-tf + {p2-tf 



— m 



{Pi ■ P2) ({t ■ e)(t ■ /i) + (t ■ e){t ■ fi)-m^{e- /x) 
-2((pi ■ t){pi ■ n){p2 ■ e) + {p2 ■ t){pi ■ e){p2 ■ /i) - {pi ■ t){p2 ■ t)(e ■ /x)) 

for unlike-helicity gluons, whereas for like-helicity gluons we have 

I^IL + = ^ ■ e)it ■ /i) + (t- e)(t ■ /i) - m2(e ■ /i)}. 

The overall factor K, is given by 



(47) 



(48) 



^'9t 9l 
3 {mtVtY 

X 



A{pi-tf + A{p2-tf -{pi-t){p2-t) 



{pi-ty{P2-ty 



b ■ V 



b ■ V 



(2e ■ z/ - m^)2 -t- (m^r^)2 (2/i ■ i/ - m^)^ ^ (m^F^)^ 



(49) 



Appendix B contains the corresponding expressions describing qq — )■ tt. 

Notice the simplicity of the matrix element squared for like-helicity gluons to top quark 
pairs. Given this simplicity and the fact that the like-helicity gluon contribution dominated 
at smaller values of the invariant mass of the tt system, it is worth exploring whether or not 
the full matrix element enhances the spin correlations in this channel. 
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The ratio of the correlated to uncorrelated^ matrix element squared, S, for like-helicity 
gluons is given by 

^_ l-^liJcorr _ mg{(t ■ e)(t ■ /i) + (t ■ e)(t ■ /i) - m,^(e ■ /i)} 

= l^i^^ / (l + /j^) + (l-/3^)Ce-^-2/32Qe-Q> 

V1 + /3V V (l-/3cte-)(l-/3cf^) 

where the last line is given in the ZMF in terms of speed of the tops, /3, and the cosine 
of the angles between t and e (cjg), i and /i (cf^) and e and /i (cg^). The range of S is 
between (2,0). At threshold, /3 — )• 0, the maximum of S occurs when the charged leptons are 
parallel, Cgfi = +1, whereas the minimum occurs when the charged leptons are back-to-back, 
Cgfj. = — 1, independent of their correlation with the top-antitop axis. 

For non-zero /3, the maximum (minimum) still occurs when the charged leptons are 
parallel (back-to-back), but they are now correlated with the top-antitop axis. The fact 
that the charged leptons are more likely to have their momenta being parallel rather than 
back-to-back is what is expected for top quark pairs that have spins which are anti-aligned, 
i.e. LL or RR. However, here the enhancement is even stronger than what one would naively 
expect because the interference between LL and RR strengthens the correlation between the 
momenta of the two charged leptons. This argument suggests looking at the AR, Ar] and 
distributions of the two charged leptons with a cut on the invariant mass of the top-antitop 
system. 



VI. CORRELATION-SENSITIVE ANGULAR DISTRIBUTIONS 
A. A(j) Distribution in Dilepton Events 

In Fig. |H]we have plotted the A0 distribution in the dilepton channel for tt production 
incorporating a cut which restricts the true invariant mass of the tt pair to less than 400 GeV. 
This plot shows results for both the fully-correlated and the uncorrelated matrix elements 
including both gg — )■ ti and qq — )■ ti channels. A clear distinction between the correlated and 
uncorrelated decays^ is seen in this figure: the difference between the two A(j) distributions 
is about 40% at both Acp = (enhancement) and Acp = tt (suppression). With this cut, 
10% of the total cross section for tt production survives at leading order. 

Unfortunately, the presence of the two neutrinos in the final state of dilepton events 
complicates the selection of events. The available kinematic constraints leave an up to 8- 
fold ambiguity^ in the reconstruction of the neutrino momenta from the available observed 



We call the decay of a top or anti-top quark into a M^-boson and 6-quark uncorrelated if this decay 
is spherical in the top quark rest frame and thus independent of the top quark spin. The M^-boson is 
then assumed to decay in the usual (fully correlated) manner. The uncorrelated matrix elements squared 
are then given by {\A\j^i^ + \A\\^)uncoTr ^ JC {t ■ e){i ■ fi) {m'^{t ■ t)} and (1^1^^ + \A\'ljj)u,icorr = 
IC{t ■ e)(? • /i){[2(pi • t)ip2 ■ t) - m2(pi •p2)][(pi • tf + (p2 • tf + ml{pi ■ P2)] I {pi ■ P2f} ■ 
^ The corresponding A?7 distribution shows almost no difference between correlated and uncorrelated matrix 

elements. Thus all the difference in the Ai? distributions comes from the A(/i distributions. 
^ The presence of a pair of quadratic constraints in the kinematic equations leads to up to 4 solutions for 
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momenta and energies. Given the ease of measuring the azimuthal angles of charged leptons, 
it is worthwhile to investigate an alternative to the true tt invariant mass cut. 

Dilepton Events 



0.5 

< 

T3 0.4 
\ 
b 

^ 0.3 
0.2 



0.1 

12 3 

FIG. 8: The differential distribution of A(p, [I/ctt) da/d{A(j)). The solid curve is for the fully 
correlated case whereas the dashed curve assumes that the top quarks decay spherically in their 
respective rest frames. A cut restricting the invariant mass of the tt pairs to a maximum of 400 GeV 
has been applied to these distributions. 

The simplest option one could imagine is to simply take the (naive) unweighted average 
{rriu) of all of the real solutions returned by the neutrino reconstruction algorithm. In Fig. [9] 
we present the results of implementing just that option: the cut used to generate this figure 
requires that {rritt) be less than 400 GeV. With this cut approximately 5% of the total 
cross section for ti production survives at leading order. This is smaller than the fraction 
passing a 400 GeV cut on the true value of mu since only those events where all the spurious 
solutions are sufficiently small will survive. On the other hand, the sample passing this cut 
will contain a few events where the true value of rriu is above 400 GeV, but, because the 
spurious solutions produced smaller values, the average was below 400 GeV. Turing to the 
A(f) distribution and comparing to the cut on the true value of rriu, one sees a rather large 
effect on the shape of the distributions. However, this effect (an enhancement near A(j) = 
and a depletion near A(f) = vr) occurs for both the correlated and uncorrelated data sets. 
Thus the difference between the two distributions remains at roughly the 40% level. No 
effort has been to optimize this invariant mass cut. Perhaps there are other variables that 
will do better than unweighted average {rriu), or perhaps 400 GeV is not the optimal cut 



any given pairing of the b jets with the two charged ieptons. Since there are two possibie pairings, as 
many as 8 different soiutions couid resuit. However, not aii of these solutions need be real, and so there 
are often fewer than the maximum possible number of solutions. 



T I I I I =1 I I I I I I I I I 

Vi" = 14 TeV : 
true m,T < 400 GeV - 



correlated 




gg + qq 

J I I I I I I I I I I I I L 
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FIG. 9: The differential distribution of Ac/), (1/(Tt) da/d{A4>). Tlie solid curve is for the fully 
correlated case whereas the dashed curve assumes that the top quarks decay spherically in their 
respective rest frames. A cut restricting the average reconstructed invariant mass of the ti pairs 
to a maximum of 400 GeV has been applied to these distributions. 

value. Nevertheless, what we have here is a proof-in-principle that these correlations can be 
measured in an experiment. 

B. ZMF cosO Distribution for Lepton-plus-jets Events 

Turning to the lepton-plus-jets channel, we have found that the cosine of the opening 
angle between the charged lepton and the d-quark jet as viewed in the zero momentum 
frame (ZMF) is sensitive to the presence or absence of correlations between production 
and decay [see also the discussion near Eq. flSOp ]. For this type of event the kinematic 
constraints provide more equations than unknowns. Thus, the ZMF may be reconstructed 
without ambiguity more than 98% of the time by discarding those solutions which do not 
pass some rudimentary quality-control cuts: the neutrino energy ought to be positive in a 
correctly reconstructed event; futhermore, the neutrino and top quark mass-shell constraints 
ought to be satisfied to sufficient accuracy.^ Discarding the ~2% of events that have more 
than one viable reconstruction of the ZMF is an acceptible option. 

For the purposes of generating this distribution, we define the ci-quark jet to be the 
jet which is spatially the closest to the 6-tagged jet in the W rest frame, as was used in 



^ Because of jet energy measurement uncertainties, it is not expected that a correctly reconstructed top 
quark or neutrino will be precisely on mass shell. The exact definition of "sufficient precision" is therefore 
a detector-dependent issue to be determined by the experimental collaborations. 
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Ref. j2|. This is equivalent to using the lowest energy jet in the top quark rest frame as 
advocated in Ref. [l5|. Fig. [10] displays the results for this distribution using only those 
events that pass the cut rriu < 400 GeV. For fully-correlated top decays this distribution is 
nearly flat, whereas for spherical decays there is a strong peaking near cos^ = —1. That 
is, because of the correlations between production and decay, the lepton and the d-jet tend 
to be significantly less back-to-back in the ZMF than if no such correlations were present. 
Approximately 9% of the total cross section for ti production survives this cut at leading 
order, even at reduced center-of-mass energy. 
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FIG. 10: The differential distribution of cos^, (l/fTj') da / d{cos 9) , where 9 is the ZMF angle 
between the charged lepton and the d-quark jet (defined to be the jet which is spatially the closest 
to the 6-tagged jet in the W rest frame; this is also the jet with the lowest energy in the top quark 
rest frame). The solid curve is for the fully correlated case whereas the dashed curve assumes that 
the top quarks decay spherically in their respective rest frames. A cut restricting the invariant 
mass of the tt pairs to a maximum of 400 GeV has been applied to these distributions. 



C. Varying the Energy of the LHC 

So far we have used 14 TeV for the energy of the LHC. However, it is now clear that 
the LHC will not reach this energy for a number of years so we have investigated what 
happens for a reduced center-of-mass energy in this section. The results we describe below 
are summarized in Fig. [TTl 

The primary result of a lower center-of-mass energy is a big reduction in the tt production 
cross section because of the reduced gluon luminosity. Panel (a) in Fig. [11] tracks the 
leading order cross section from 7 to 14 TeV center-of-mass energy. We see that a factor 
of 2 reduction in y/s produces a reduction of about a factor of 5 in the tt production cross 
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FIG. 11: Effects of varying the machine center-of-mass energy ^/s. (a) Total leading order cross 
section for pp — )• tt. These values should be multiplied by the branching fraction to dileptons (4.6%) 
or lepton plus jets (29%), as appropriate. We include only the e and /i channels, (b) Fraction of 
dilepton and lepton plus events with m^f < 400 GeV. For dilepton events (crosses) we employ 
the unweighted average of the up to 8 solutions for the tt invariant mass. For lepton+jets events 
(diamonds) the true value of m^f may be reconstructed and used in event selection, (c) Half of 
the area between the appropriate unit-normalized angular distributions for the fully correlated and 
spherical cases. For lepton+jets events (crosses), we use the distribtuion in cos 6, where 6 is the 
angle between the charged lepton and the d-jet candidate in the zero momentum frame of the 
event. For dilepton events (diamonds), we use the azimuthal opening angle Acj) between the two 
charged leptons. 

section. Panel (b) illustrates the fact that the fraction of dilepton and lepton plus jet events 
surviving the mu or {rriu) cut we advocate does not change very much as a/s is varied 
between 7 and 14 TeV. Finally, in panel (c) of each figure we compare the quantity A1/2, 
which is defined as half of the area between the correlated and uncorrelated predictions for 
the unit-normalized angular distributions. This quantity ranges from in the case where 
the distribution is completely independent of whether or not correlations are present up 
to 1 in the case where the two distributions (correlated decays or spherical decays) do not 
overlap at all. In both channels, there is relatively little dependence on ^/s for this measure 
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of the difference between tlie correlated and splierical cases. Thus, the biggest issue related 
to the observation of these spin correlations at the LHC running at reduced energy comes 
from the greatly diminished cross section: the correlations themselves remain at roughly the 
same level. Fortunately, even a reduction of the number number of tt pairs estimated in 
the introduction by a factor of 5 leaves ample statistics to hope for at least a preliminary 
observation of these spin correlation effects, even at reduced center-of-mass energy. 



D. NLO Effects 

Higher-order QCD effects enhance the total cross section especially near threshold. How- 
ever, previous studies on e"'"e~ — )• tt demonstrate that such corrections to the spin correla- 
tions are small, see Ref. [l6|. One can understand this physically since the emission of soft 
gluons from a top quark cannot flip the spin of the top quark. We have done a preliminary 



study of the NLO effects using both MCFM, [171] and MC@NLO, [18[, incorporating a cut 



which restricts the invariant mass of the tt pair to be less than 400 GeV. Both these Monte 
Carlos show that the tree-level effects discussed earlier in this paper are also present at 
NLO.'^ We need not check the Monte Carlo of Ref. ^ since it agrees with the other two 
Monte Carlos for the correlated decays. A detailed study at NLO where the invariant mass 
of the ti pair is reconstructed from the decay products is beyond the scope of this work. 
A NLO Monte Carlo with a switch that allows the user to toggle between fully correlated 
top quark decays and spherical decays of the tops would be very useful for such a study. 
However, no such Monte Carlo exists at present. 



VII. SUMMARY AND CONCLUSIONS 

In this paper we have shown how to observe spin correlations in top quark pair production 
at the LHC. To our surprise, the observation of these correlations is easier at the LHC than 
at the Tevatron. The reason for this is that at the LHC the dominant production mechanism 
for top quark pairs is gluon-gluon fusion, which at low s is dominated by the fusion of like 
helicity pair gluons. The fusion of like helicity gluons produces top quark pairs in a LL 
or RR helicity configuration. When such top quarks decay, they produce charged leptons 
which possess very strong azimuthal correlations. These correlations can be easily seen in 
the laboratory frame once a cut on the invariant mass of the top quark pair is made. There 
is not need to reconstruct the top quark rest frame as is required to see the correlations of 
top quark pairs at the Tevatron. The analysis has been extended to the lepton plus jets 
channel by "identifying" which jet is the ci-quark jet from the ly-boson decay on a statistical 
basis. Apart from the reduction of the total cross section, the size of these spin correlations 
is approximately independent of the energy of the LHC in the 7 to 14 TeV range. Thus, we 
expect that these effects could be observed in early running of the LHC. 



^ In Fig. 1 of Ref. [18[ one can see the size of the correlations without the invariant mass cut. 
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Appendix A: Universal Properties of the Spin Amplitudes 

In this Appendix we demonstrate the properties of the spin amplitudes and their rela- 
tionship to the Wigner d-functions and the helicity amplitudes. Of particular importance 
is the fact that spin amplitudes at either end of the spin chain contain all the information 
necessary to obtain all the other spin amplitudes by simple algebraic manipulations. This 
fact was used extensively earlier in this paper. 



1. Wigner d-functions and connection to the spin amplitudes 

The amplitude for production of a single massive particle with spin j and spin-projection 
m in the generalized spin basis described in Sect. Ill Al and illustrated in Fig. [1] may be 



written as a linear superposition of the matrix elements of the rotation operator [20 

m' 

The c?;^ appearing in Eq. (lAip are the Wigner d-functions, chosen to conform to the 
conventions of Rose 2l|. In particular, we write — ^ since our angle ^ is measured in the 
clockwise direction whereas the angle /3 in Ref. M is counterclockwise. Since 



<w(0) = (A2) 

we see that 

i^™ = A-,™(0)- (A3) 

That is, the coefficients Hm are the conventional helicity amplitudes for the process with a 
particular choice of relative phases. 

Similarly, the production of a pair of massive particles (spins ji and j2, spin projections 
mi and 1712 along independent spin axes oriented at the clockwise angles ^ and C,' with respect 
to the recoil direction in the scattering plane) may be decomposed as 

m m' 

The extension to more than two particles in the final state or to spin axes which point out 
of the production plane, involving the introduction of the Wigner D-functions in place of 
the (simpler) ci-functions, is straightforward, but beyond the scope of this appendix. 
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2. ,^ — )• ^ ± TT rule 

Intuitively, we expect that the probability for producing spin projection +m along some 
axis ought to be equal to that for producing spin projection —m along minus that axis: 

\A,,mi^±n)\ = \A,,-mm- (A5) 

That is, flipping the z component of the spin and rotation by ±7r give the same results up 
to a phase. This feature of the spin amplitudes is a consequence of the properties of the 
Wigner ci-functions appearing in Eqs. flAip and flA4p . 



3. Differential relations among the spin amplitudes 

Consider a state with total spin j and projection jz = m in the general spin basis 
illustrated in Fig. 1. Applying the rotation operator converts this to a state where the spin 
axis is at ,^ + A,^ instead |20(] : 

Aj,U^ + AO = e^^«^M,,^(0. (A6) 
In the limit A,^ — t- Eq. flA6p may be rewritten as 

^^jAO = liJ+ - J-)A,AO- (A7) 

In this expression, we have replaced Jy by the appropriate linear combination of raising and 
lowering operators. Thus 

d 1 

- ^y/j{j + l)-m{m-l)A,,m-iiO (A8) 



4. Starting at top or bottom of spin chain 



It is useful to record the explicit results of applying Eq. (lASP to the ends of the spin 
chain: 



A±(i-i)(0 = T 

A second differentiation allows us to conclude that 

1 



A 



■i,±0-2) 



(0 



(A9) 



(AlO) 



This process could be repeated as many times as required. However, when used in conjunc- 
tion with the ^ — )■ ^ ± vr rule of Eq. ( lASp . Eqs. ( 1A9P and ( lAlOp allow the calculation of all 
of the amplitudes in the spin chain for spins up to and including |, starting from either end 
of the chain (m = ±j). The relative simplicity of the operations involved in these relations 
makes for a substantial computational savings over the calculation of the complete set of 
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amplitudes by direct means. Of course, once you have all the spin amplitudes then the 
helicity amplitudes can be easily obtained, including relative phases, by setting ^ = 0. 

Thus, the spin amplitudes at the top and bottom of the spin chain contain all the in- 
formation about a given process and all other spin amplitudes and helicity amplitudes can 
be derived from them. This special property of these spin amplitudes is simply reflected 
in the Wigner (i-functions, d-j^ for m = {—j,---,j}, which form a set of (2j+l) linearly 



independent functions [22 



5. Examples: 

An explicit illustration of the relationships contained in Eqs. (1A9P and (lAlOp is provided 
by the gg — > ti process considered in this paper. After setting ^' = ^ (back-to-back spin 
axes in the ZMF), we can organize the amplitudes according to the total spin in the final 
state. In doing this we need to keep in mind that for this choice of spin axes, unlike spin ti 
pairs have spins that point in the same spatial direction. Thus 



and 



IAi(OI = lAigg^t^k 



|^i,-i(OI = lAigg^^i^ 



(All) 



(A12) 



The ^-dependent hnear combination of the and H amplitudes must be Aifi, the orthog- 
onal (^-independent) combination is ^o,o: 



|A,o(OI = ^\^(99 ^ h^) + Aigg ^ t^k)l 
|A,o(OI = ^ ^t^t) - A{gg ^ tj;)|. 



(A13) 
(AM) 



For like-helicity gluons Eqs. (31) and (32) lead to 



|A,o(OI ~ v^7"Vcd 
IA-i(OI ~ T^V^d 

The three j = 1 amplitudes in (1A15I) satisfy 

IA,o(OI = 

and 



and |A,o(OI~v^' 



'7 



-1 



1 + 2 



AiAO 



(A15) 



(A16) 



(A17) 



as implied by Eqs. (^M^ and (IXTOj) . 

A different realization of these relationships is provided by the unlike-helicity gluons for 
which Eqs. (22) and (23) lead to 



|^i,i(OI ~ /3se(l + cecj + 7 ^S0S^) 
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|A,o(OI ~ V2/3se\ceSi: 



-1, 



7 sec^l 



and |^o,o(OI = 0. 



(A18) 



This is identical (up to overall factors) to what happens for the processes qq — )■ tt and similar 
to what happens for e^e^ — )■ ti [5]. The spin amplitudes for e"^e~ — )■ Z/i 23|, |2J| also satisfy 
Eqs. (IXTej) and (IXTTl) . 

Finally, we have verified that the processes e'^e^ — > W^W~ and provide examples 
of the J = 2 versions of Eqs. (lA9p and (lAlOp . Indeed, the derivative relations between the 
spin amplitudes for these processes noted in Ref. [23| are a direct consequence of Eq. (]A9ll . 



Appendix B: The process qq — > M 

For completeness we give here the matrix element squared for qq — )■ ti with the subsequent 
decay of the top quarks. Starting from Eqn. ^ and using the substitutions given in Eq. (jlHD, 
it is easy to add the decays of the on-mass-shell top quarks {t — )■ h + e + v and t — >• b + ix + u). 
Thus, the total matrix element squared for top quark production and decay via quark- 
antiquark annihilation, summed over the colors of the incoming and outgoing quarks, is 
given by 



I^IL + I^I 



2 

LR 



/C,5<|2(t-e)(t-/i) [{q-tf + {q-t) 

{q-q)({t-e){t- fi) + {i ■ e) {i ■ fx) - {e ■ fi)^ 
-2 f (g ■ t) {q ■ /x) (g ■ e) + (g ■ t) {q ■ e) (g ■ /x) - (g ■ t) (g ■ t){e ■ /i) 



-m. 



(Bl) 

This has the same functional form as the part of Eqn. (H7|) in the second set of curly brackets. 
The overall factor ICqq is given by 



6„4 



9l 



h-u 



b ■ V 



(g ■ g)2 {m^TtY (2e ■ z/ - mlf + {m^T^f (2/x ■ z7 - mlf + {m^T^f ' 
Eqs. f lBip and f lB2p are the Lorentz-invariant equivalents of Eqs. (4) and (5) of Ref. 



(B2) 
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See talk by Florencia Canelli at Lepton Photon 2009, http://lp09.desy.de/ 
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